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Abstract 

A new set of measurements of the top quark mass are presented, based on 
the proton-proton data recorded by the CMS experiment at the LHC at - 
8 TeV corresponding to a luminosity of 19.7 fb _1 . The top quark mass is mea¬ 
sured using the lepton+jets, all-jets and dilepton decay channels, giving values 
of 172.35 ± 0.16 (stat) ± 0.48 (syst) GeV, 172.32 ± 0.25 (stat) ± 0.59 (syst)GeV, and 
172.82 ± 0.19 (stat) ± 1.22 (syst) GeV, respectively. When combined with the pub¬ 
lished CMS results at ^/s = 7 TeV, they provide a top quark mass measurement of 
172.44 ± 0.13 (stat) ± 0.47 (syst) GeV. The top quark mass is also studied as a function 
of the event kinematical properties in the lepton+jets decay channel. No indications 
of a kinematic bias are observed and the collision data are consistent with a range of 
predictions from current theoretical models of tt production. 


Published in Physical Review D as doi : 10.1103/PhysRevD.93.072004. 


© 2016 CERN for the benefit of the CMS Collaboration. CC-BY-3.0 license 


*See Appendix[A]for the list of collaboration members 









1 


1 Introduction 


The mass of the top quark (m t ) is one of the fundamental parameters of the standard model 
(SM). A precise measurement of its value provides a key input to global electroweak fits and to 
tests of the internal consistency of the SM [U0- Its value leads to constraints on the stability of 
the electroweak vacuum 00 and affects models with broader cosmological implications 0111. 

The most precise measurements of m t have been derived from combinations of the results from 
the CDF and DO experiments at the Tevatron, and ATLAS and CMS at the CERN LHC. The 
current combination from the four experiments gives a top quark mass of 173.34 ± 0.76 GeV |7j, 
while the latest combination from the Tevatron experiments gives a mass of 174.34 ± 0.64 GeV fO 
The Tevatron combination is currently the most precise measurement and it includes all of the 
current Tevatron measurements. In contrast, the current four experiment combination has not 
been updated since 2013 and does not include the latest Tevatron and LHC measurements, in 
particular the measurement from ATLAS using a combination of the lepton+jets and dilepton 
channels H9j. 

Beyond the leading order (LO) in quantum chromodynamics (QCD), the numerical value of 
depends on the renormalization scheme llOl flT 'l. The available Monte Carlo (MC) generators 
contain matrix elements at LO or next-to-leading order (NLO), while higher orders are approx¬ 
imated by applying parton showering. Each of the measurements used in the combinations 
has been calibrated against the mass implemented in a MC program. Given the precision of the 
experimental results, a detailed understanding of the relationship between the measurements 
and the value of m t in different theoretical schemes is needed. Current indications are that the 
present measurements based on the kinematic reconstruction of the top quark mass correspond 
approximately to the pole ("on-shell") mass to within a precision of about 1 GeV liT2l . 

At the LHC, top quarks are predominantly produced in quark-antiquark pairs (tt) and top 
quark events are characterized by the decays of the daughter W bosons. This leads to exper¬ 
imental signatures with two jets associated with the hadronization of the bottom quarks and 
either a single lepton (e, }i), one undetected neutrino and two light quark jets (lepton+jets chan¬ 
nel), or four light quark jets (all-jets channel), or two leptons (ee, e}i, }iu) and two undetected 
neutrinos (dilepton channel). While the events which contain leptonic r decays are included 
in the analysis samples, they contribute very little to the mass measurements as their yields 
are negligible. The results presented in this paper focus on the analysis of data in these three 
channels recorded by the CMS experiment in the 2012 part of what is commonly referred to as 
Run 1 of the LHC. 


The paper is organized as follows. The main features of the detector and the data are dis¬ 
cussed in Sections [2] and [3] Section [4] is a discussion of the analysis techniques, which lead to 
the measurements of Section [5] The categorization of the systematic uncertainties is presented 
in Section [6j followed by the full results for the three decay channels in Section [7J Section [8] 
presents a study of m t as a function of the kinematical properties of the tt system in the lep¬ 
ton+jets channel. This is followed in Sections [9] and 10 which discuss the combination of the 
measurements and the final result for m t , respectively. 


2 The CMS detector 

The central feature of the CMS apparatus is a superconducting solenoid of 6 m internal diam¬ 
eter, providing a magnetic field of 3.8 T. Within the solenoid volume are a silicon pixel and 
strip tracker, a lead tungstate crystal electromagnetic calorimeter (ECAL), and a brass and scin- 
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3 Data sets 


tillator hadron calorimeter (HCAL), each composed of a barrel and two endcap sections. The 
tracker has a track-finding efficiency of more than 99% for muons with transverse momentum 
pj > lGeV and pseudorapidity \tj\ < 2.5. The ECAL is a fine-grained hermetic calorimeter 
with quasi-projective geometry, and is distributed in the barrel region of \t]\ < 1.48 and in 
two endcaps that extend up to \rj\ < 3.0. The HCAL barrel and endcaps similarly cover the 
region ;/| < 3.0. In addition to the barrel and endcap detectors, CMS has extensive forward 
calorimetry. Muons are measured in gas-ionization detectors, which are embedded in the steel 
flux-return yoke outside of the solenoid. A more detailed description of the CMS detector, to¬ 
gether with a definition of the coordinate system used, and the relevant kinematic variables 
can be found in Ref. Ifl3| . 

3 Data sets 

The measurements presented in this paper are based on the data recorded at a center-of-mass 
energy of 8 TeV during 2012, and correspond to an integrated luminosity of 19.7 fb -1 . 

3.1 Event simulation and reconstruction 

Simulated tt signal events are generated with the MadGraph 5.1.5.11 LO matrix element gen¬ 
erator with up to three additional partons Iil4l . MadSpin ffl5l is used for the decay of heavy 
resonances, PYTHIA 6.426 for parton showering jl6l using the Z2* tune, and TAUOLA fl7l for de¬ 
cays of r leptons. The most recent PYTHIA Z2* tune is derived from the Z1 tune lfl8l . which uses 
the CTEQ5L parton distribution function (PDF) set, whereas Z2* uses CTEQ6L fl9l . A full sim¬ 
ulation of the CMS detector based on Geant4 11208 is used. The tt signal events are generated 
for seven different values of m t ranging from 166.5 to 178.5 GeV. The W/Z+jets background 
events are generated with MadGraph 5.1.3.30. The diboson background (WW, WZ, ZZ) is 
simulated using PYTHIA 6.426 using the Z2* tune. The single top quark background is simu¬ 
lated using POWHEG 1.380 ff2lT(25ll , assuming an m t of 172.5 GeV. The tt, W/Z+jets, and single 
top quark samples are normalized to the theoretical predictions described in Refs. M26H30fl . The 
simulation includes the effects of additional proton-proton collisions (pileup) by overlapping 
minimum bias events with the same multiplicity distribution and location as in data. 

Events are reconstructed using a particle-flow (PF) algorithm j|3Tl 321. This proceeds by recon¬ 
structing and identifying each final-state particle using an optimized combination of all of the 
subdetector information. Each event is required to have at least one reconstructed collision 
vertex. The primary vertex is chosen as the vertex with the largest value of £ Pj of the tracks 
associated with that vertex. Additional criteria are applied to each event to reject events with 
features consistent with detector noise and beam-gas interactions. 

The energy of electrons is determined from a combination of the track momentum at the 
primary vertex, the corresponding ECAL energy cluster, and the sum of the reconstructed 
bremsstrahlung photons associated with the track ff33l . The momentum of muons is obtained 
from the track momentum determined in a combined fit to information from the silicon trackers 
and the muon detectors Il34l . The energy of charged hadrons is determined from a combination 
of the track momentum and the corresponding ECAL and HCAL energies, corrected for the 
suppression of small signals and calibrated for the nonlinear response of the calorimeters. Fi¬ 
nally, the energy of the neutral hadrons is obtained from remaining calibrated HCAL and ECAL 
energies. As the charged leptons originating from top quark decays are typically isolated from 
other particles, a relative isolation variable (T re l) is constructed to select lepton candidates. This 
is defined as the scalar sum of the pj values of the additional particles reconstructed within 
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an angle A R = \J (A?/) 2 + (A (p) 2 of the lepton direction, divided by the pj of the lepton. Here 
A?/ and A (p are the differences in the pseudorapidity and azimuthal angles between the lepton 
direction and other tracks and energy depositions. A muon candidate is rejected if I re \ > 0.12 
for A R = 0.4, and an electron candidate is rejected if I rc \ > 0.10 for A R = 0.3. 

Jets are clustered from the reconstructed PF candidates using the anti-kj algorithm Il35l with 
a distance parameter of 0.5, as implemented in the FASTjET package Ii36l . The jet momen¬ 
tum is determined from the vector sum of the momenta of the particles in each jet, and is 
found in simulation to be within 5 to 10% of the jet momentum at hadron level for the full 
pj range If37l . Corrections to the jet energy scale (JES) and the jet energy resolution (JER) are 
obtained from the simulation and through in-situ measurements of the energy balance of ex¬ 
clusive dijet, photon+jet, and Z+jet events. Muons, electrons, and charged hadrons originating 
from pileup interactions are not included in the jet reconstruction. Missing transverse energy 
(Ej' ss ) is defined as the magnitude of the negative vector pj sum of all selected PF candidates 
in the event. Charged hadrons originating from pileup interactions are not included in the cal¬ 
culation of £j lss . Jets are classified as b jets through their probability of originating from the 
hadronization of bottom quarks, using the combined secondary vertex (CSV) b tagging algo¬ 
rithm, which combines information from the significance of the track impact parameters, the 
kinematical properties of the jets, and the presence of tracks that form vertices within the jet. 
Three different minimum thresholds are used for the CSV discriminator to define the loose 
(CSVL), medium (CSVM), and tight (CSVT) working points. These have b tagging efficiencies 
of approximately 85%, 67%, and 50%, and misidentification probabilities for light-parton jets 
of 10%, 1%, and 0.1%, respectively l38l . 

3.2 Event selection 

For the lepton+jets channel we use the data collected using a single-muon or single-electron 
trigger with a minimum trigger pj threshold for an isolated muon (electron) of 24 GeV (27 GeV), 
corresponding to an integrated luminosity of 19.7 fb 1 . We then select events that have exactly 
one isolated muon or electron, with pj > 33GeV and |//| < 2.1. In addition, at least four jets 
with pj > 30 GeV and | ?/1 < 2.4 are required. Jets originating from b quarks (denoted as b jets) 
are identified using the CSV algorithm at the medium working point 15SI . With the require¬ 
ment of exactly two b-tagged jets among the four jets with the highest pj, 104 746 tt candidate 
events are selected in data. From simulation, the sample composition is expected to be 93% 
tt, 4% single top quark, 2% W+jets, and 1% other processes. Figure [l] shows the comparison 
of the data and simulation for the selected events in some representative distributions. The 
simulation shown is not corrected for the uncertainty in the shape of the top quark pj dis¬ 
tribution 1391 . which accounts for almost all of the slope visible in the data/MC ratio plots. 
However, even without making a correction, the data and simulation are consistent within the 
quoted uncertainties. 

For the all-jets channel we use the data collected using a multijet trigger, corresponding to 
an integrated luminosity of 18.2 fb '. The trigger requires the presence of at least four jets, 
reconstructed from the energies deposited in the calorimeters, with transverse momenta pj > 
50 GeV. Since fully hadronically decaying top quark pairs lead to six partons in the final state, 
events are required to have at least four jets with pj > 60 GeV and a fifth and sixth jet with 
pj > 30 GeV. Jets originating from b quarks are identified using the CSV b algorithm at the 
tight working point Il38l . With the requirement of exactly two b-tagged jets among the six 
leading ones, 356 231 candidate events are selected. From simulation, the sample is expected 
to be dominated by the QCD multijet background and to have a signal fraction of about 13%. 
The QCD multijet background cannot be reliably simulated and we determine its kinematic 
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Figure 1: Distributions for the lepton+jets channel of (upper left) lepton pj, (upper right) 
missing transverse energy, (lower left) leading jet pj, (lower right) second-leading jet pj for 
data and simulation, summed over all channels and normalized by luminosity The vertical 
bars show the statistical uncertainty and the hatched bands show the statistical and systematic 
uncertainties added in quadrature. The lower portion of each panel shows the ratio of the 
yields between the collision data and the simulation. 
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dependence from a control sample in the data. The background normalization is determined 
as a part of the fit process, which is discussed in Section [T2] 

For the dilepton channel, events are required to pass the triggers appropriate for each of the 
three channels. The ep channel uses a logical OR of two triggers that require a muon of pj > 
17 (or 8) GeV and an isolated electron of 8 (or 17) GeV. Dimuon events must pass a trigger 
which requires a pj > 17 GeV for the muon with the highest ("leading") pj and 8 GeV for the 
second-leading muon. Similarly, dielectron events must satisfy a trigger with a threshold of 
pj > 17 GeV for the leading electron and 8 GeV for the second-leading electron. In this case 
both electrons are required to be isolated. In all three cases the amount data corresponds to a 
luminosity of 19.7ftC 1 . We select events for analysis if they have two isolated opposite-sign 
leptons with pj > 20 GeV and |f/| < 2.4 (2.5) for muons (electrons). Jets originating from 
b quarks are identified with the CSV algorithm at the loose working point 1381 . Events are 
retained if they have at least two b-tagged jets. Background contamination from low-mass 
resonances is reduced by demanding a dilepton pair invariant mass, niff, of at least 10 GeV. To 
suppress the background from Z boson decays, events with ee and pp signatures are required 
to have E™ ISS > 40 GeV, and to fall outside of the dilepton invariant mass window 76 < niff < 
106 GeV. The remaining Drell-Yan background is estimated from the data using the ratio of the 
event yield inside vs. outside the invariant mass window Il38l . After all of the requirements, we 
find 41125 candidate events in data for which the sample compositions is expected to be 95% tt, 
3% single top quark, 2% Drell-Yan, and <0.3% other processes. Figure [2]shows a comparison of 
the data and simulation for events with at least one b jet for some representative distributions. 
As with the lepton+jets plot (Fig. |T|) the simulation is not corrected for the discrepancy in the 
top quark pj distribution, leading to the slopes visible in the data/MC plots. 


4 Analysis techniques 

The measurements discussed in the following sections use analysis techniques in which either 
m t alone is determined or m t and the overall jet energy scale factor (JSF) are determined simul¬ 
taneously. For the lepton+jets and the all-jets channels we use analyses based on the ideogram 
technique (Section [4/1) . While the ideogram technique provides the most precise measure¬ 
ments, it is not suitable for dilepton events where the presence of more than one neutrino 
introduces uncertainties in the use of the measured E™ lss . Instead, for the dilepton channel, we 
use the Analytical Matrix Weighting Technique (AMWT) method (Section |4.2[ ). 

4.1 One- and two-dimensional ideogram analyses 

The ideogram method is a joint maximum likelihood fit that determines m t and, optionally, the 
JSF from a sample of selected tt candidate events in the lepton+jets or all-jets channels. The 
observable used for measuring m t is the mass m [ ll estimated by a kinematic fit [[401 . The kine¬ 
matic fit constrains the candidates for the tt decay products to the hypothesis of the production 
of two heavy particles of equal masses, each one decaying to a W boson and a bottom quark, 
where the W boson invariant mass is constrained to 80.4 GeV I14T1 . The JSF is defined as a mul¬ 
tiplicative factor to be applied in addition to the standard jet energy corrections (JEC) l37l to 
the four-momenta of the jets. The JSF is determined from the invariant masses of the jet pairs, 
m(“°, associated with the W bosons before the jet momenta are constrained by the kinematic 
fit. For the case of a simultaneous fit to both m t and the JSF (2D approach), no prior knowledge 
of the JSF is assumed. If only m t is fitted (ID approach), the jet energy scale determined from 
the JEC is taken as the JSF prior, fixing it to unity. A third category of fits (hybrid approach) 
incorporates the prior knowledge about the jet energy scale by using a Gaussian constraint. 
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Figure 2: Distributions for the dilepton channel: (upper left) leading lepton pj, (upper right) 
second-leading lepton pj, (lower left) leading jet pj, (lower right) second-leading jet pj for 
data and simulation, summed over all channels and normalized by luminosity. The vertical 
bars show the statistical uncertainty and the hatched bands show the statistical and systematic 
uncertainties added in quadrature. The lower portion of each panel shows the ratio of the 
yields between the collision data and the simulation. 
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P(JSF), centered at 1 with a variance depending on the total JEC uncertainty. For the hybrid 
analysis in the lepton+jets channel, the JSF determined from the Wboson decays and the jet 
energy scale from the JEC are given equal weight in the fit. In contrast, for the hybrid fit in the 
all-jets channel, the jet energy scale from the JEC contributes 80% of the information, because 
of the larger uncertainty on the JSF from the 2D fit. 

The distributions of m[ d and m((', co are obtained from simulation for three to seven different m t 
and three to five different JSF values for the tt signal, and from simulated background events 
(lepton+jets) or the control sample for the multijet background (all-jets). From these distri¬ 
butions, probability density functions are derived separately for different cases of jet-parton 
assignments for the signal, and for the background contribution. The signal functions depend 
on m t and JSF, and are labeled P(7n^ t \ni t ,]S¥) and P(m^. C0 |nq,JSF), respectively, for an event in 
the final likelihood. 

The likelihood for measuring nit and the JSF in an observed data sample can be expressed as: 

£ (sample|mt, JSF) = C (event|m t , JSF) li ’ event , (1) 

events 

where the event weight m event = c Ya=\ Pgof (0 is used in the lepton+jets analysis to reduce 
the impact of events for which the chosen permutation of the jets is incorrect. Here, c is a 
normalization constant and the remaining quantities are defined as in Eq. (|2jl. For the all-jets 
channel, m ev ent = 1 is used. The event likelihoods (or ideograms) are given by 

C (event | m t , JSF) = J2 P s° f (0 j/sigh’sig K,JSf) + (l -fag) Pbkg (mg,m^ > )|, 

( 2 ) 

where the index i runs over the n selected permutations of an event that each have a goodness- 
of-fit probability P go f assigned from the kinematic fit. The signal fraction f s j g is assumed to be 1 
for the lepton+jets channel and is left as a free parameter of the fit for the all-jets channel. The 
background term g is independent of both m t and the JSF for backgrounds determined from 
the collision data. 

As the W boson mass is fixed to 80.4 GeV in the fit, the observables m[' 1 and m(y C0 have a low 
correlation coefficient (less than 5%) and the probability density P can be factorized into one¬ 
dimensional expressions, 

P [mf, mJy co | m u JSf) = Yj fi P i I m n J SF ) p j ( m w°I"G JSF), (3) 

i 

where the index j denotes the different jet-parton permutation classes defined for the measure¬ 
ment. Their relative fraction fj is either determined from the simulated sample with m t/gen = 
172.5 GeV or by the fit. 

The most likely m t and JSF values are obtained by minimizing —2 In C (sample|m t , JSF) for the 
2D and hybrid analyses. For the ID analyses only nit is determined and the JSF is set to unity 
during the minimization. 

4.2 Analytical matrix weighting technique (AMWT) 

The measurement of m t for the dileptonic tt decays is performed using the analytical matrix 
weighting technique (AMWT). This is based on a matrix weighting technique used by the DO 
collaboration Il42l . combined with an analytical algorithm to find solutions of the kinematic 



8 


4 Analysis techniques 


equations Il43l . The method allows the determination of m t with the assumption of JSF = 1, 
and in this sense, the results are comparable to the ID fits performed in either the lepton+jets 
or all-jets channels (see Section [44) . 

In dileptonic tt decays, the final state consists of two charged leptons, two neutrinos, and two 
b quarks, resulting in 18 unknowns: three momentum components for each of the six final 
state particles. Of these, we observe the momenta of the two charged leptons, the momenta of 
the two jets, and the momenta of all of the other charged particles and jets. If there are more 
than two jets in an event we have to select the jets to assign to the b quarks from the decay of 
the top quark pair. We preferentially assign b-tagged jets to these. Hence, after physics object 
reconstruction, we measure the following observables for each event: 

• the momenta p£+ and of the charged leptons from the W + and W decays, 

• the momenta and of the b and b quarks produced by the t and t quark decays, 

• the total transverse momentum pj of the tt pair. 

This leaves four unknowns that must be solved analytically. Conservation of four-momentum 
provides the following four constraints on the kinematics, if a hypothetical value for the top- 
quark mass is assumed: 

• the masses trif i v and nif y of the lepton-neutrino pairs from the W + and W decays 
are constrained to be 80.4 GeV ll4H . 

• the masses of the systems of particles from the t and t decays must equal the hypoth¬ 
esized mass of the top quark. 

Hence, the system of equations is appropriately constrained. However, there is not a unique so¬ 
lution, because the equations are nonlinear. For a given assignment of reconstructed momenta 
to final-state particles there can be up to four solutions for the neutrino momenta such that the 
event satisfies all of the constraints. There is a twofold ambiguity of assigning jet momenta to 
the b and b jets, which doubles this to eight possible solutions. We follow the algorithm given 
in Refs. [44:, 451 to find these solutions. In rare cases, a latent singularity in the equations used 
to find these solutions can prohibit the calculation of the longitudinal momenta. In such events, 
a numerical method is employed to find the incalculable variables Ii44l . 

For each event, we find all solutions of neutrino momenta for hypothesized top quark masses 
between 100 and 600 GeV, in 1 GeV increments. In general, we expect solutions to be found for 
a large range of mass hypotheses. To each solution we assign a weight w given by 06): 


w(X\m t ) = 


E f (*l) f (*2) 

initial partons 


p (E e +\m t ) p {E e -\m t ), 


(4) 


where X represents the momentum vectors of the final state particles as obtained from the 
solutions of the kinematic equations. We sum the parton distribution functions F(x), evaluated 
at Q 2 = m 2 , over the possible LO initial parton states (uu, uu, dd, dd, and gg); X] and X 2 are 
the Bjorken x values for the initial-state partons which can be computed from the momenta of 
the final-state particles. The function p(E\m t ) is the probability density of observing a charged 
lepton of energy E in the rest frame of a top quark of mass m t , given by |46): 

Am t E(m 2 — — 2m t E) 

PE '” h) = + (5) 

where the b quark mass, ?»b, is set to 4.8 GeV, and the W boson mass, Mw, to 80.4 GeV. For 
each »z t hypothesis, we find an overall weight by summing the weights of all solutions found. 
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To compensate for mismeasurements of the momenta due to finite detector resolution or the 
loss of correlation between the jet and quark momentum because of hard-gluon radiation, we 
account for the jet energy resolution during reconstruction. Every event in both the collision 
and simulated data is reconstructed 500 times, each time with jet momenta drawn randomly 
from a Gaussian distributions of widths given by the detector resolution and with means given 
by the measured momenta. After this randomization procedure, approximately 96% of all 
events in both the collision and simulated data have at least one solution, and hence a top 
quark mass estimator. The final weight curve of each event is given by the average of the 
weight distributions from each of the 500 randomizations, after excluding the cases for which 
there is no valid solution. This distribution serves as a measure of the relative probability that 
the observed event occurs for any given value of m t . 

The estimator for m t is then the hypothesized mass with the highest average sum weight for 
each event, called the AMWT mass, /n^ MWT . 


5 Mass measurements 

5.1 The lepton+jets channel 

To check the compatibility of an event with the tt hypothesis and improve the resolution of 
the reconstructed quantities, a kinematic fit |j40l is applied to the events. For each event, the 
inputs to the fitter are the four-momenta of the lepton and the four leading jets, the missing 
transverse energy, and their respective resolutions. The fit constrains these to the hypothesis of 
the production of two heavy particles of equal mass, each one decaying to a Wboson with an 
invariant mass of 80.4 GeV BTI and a bottom quark. It minimizes x 2 = (x — x m ) T E 1 (x — x m ) 
where x m is the vector of measured observables, x is the vector of fitted observables, and E 1 is 
the inverse error matrix which is given by the resolutions of the observables. The two b-tagged 
jets are candidates for the bottom quarks in the tt hypothesis, while the two untagged jets serve 
as candidates for the light quarks for one of the Wboson decays. This leads to two possible 
parton-jet assignments per event and two solutions for the z component of the neutrino mo¬ 
mentum. 

For simulated tt events, the parton-jet assignments are classified as correct permutations, wrong 
permutations, and unmatched permutations. The correct permutation class includes those events 
for which all of the quarks from the tt decay (after initial-state parton shower) are correctly 
matched to the selected jets within a distance A R < 0.3. The wrong permutations class cov¬ 
ers the events for which the jets from the tt decay are correctly matched to the selected jets, 
but where two or more of the jets are interchanged. Lastly, the unmatched permutations class 
includes the events for which at least one quark from the tt decay is not matched unambigu¬ 
ously to any of the four selected jets. To increase the fraction of correct permutations, we re¬ 
quire Pg 0 f > 0.2 for the kinematic fit with two degrees of freedom. This selects 28 295 events 
for the mass measurement, with an estimated composition of 96.3% tt signal and 3.7% non-tt 
background, which is dominated by single top quark events. In the mass extraction, the per¬ 
mutations are weighted by their Pgoi values, and the effective fraction of correct permutations 
among the tt signal improves from 13% to 44%, while the fractions of wrong and unmatched 
permutations change from 16% to 21% and 71% to 35%, respectively, determined in simulation. 

Figure |3] shows the distributions before and after the kinematic fit and ^gof selection of the 
reconstructed mass m^ co of the Wboson decaying to a qq pair and the mass m[ eco of the corre¬ 
sponding top quark for all possible permutations. 
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Figure 3: Reconstructed masses of (upper left) the W bosons decaying to qq pairs and (upper 
right) the corresponding top quarks, prior to the kinematic fitting to the tt hypothesis. Panels 
(lower left) and (lower right) show, respectively, the reconstructed W boson masses and the 
fitted top quark masses after the goodness-of-fit selection. The total number of permutations 
found in simulation is normalized to be the same as the total number of permutations observed 
in data. The vertical bars show the statistical uncertainty and the hatched bands show the 
statistical and systematic uncertainties added in quadrature. The lower portion of each panel 
shown the ratio of the yields between the collision data and the simulation. 
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CMS Lepton+jets, 19.7 fb' 1 (8 TeV) 



Figure 4: The two-dimensional likelihood (—2A log (£)) for the lepton+jets channel for the 2D, 
hybrid, and ID fits. The thick (thin) ellipses correspond to contours of —2Alog (£) = 1 (4) 
allowing the construction of the one (two) u statistical intervals of m t . For the ID fit, the thick 
and thin lines correspond to the one and two a statistical uncertainties, respectively. 

The ideogram method (Section |4.1| ) is calibrated for each combination of the top quark mass 
hypothesis, mf en and JSF values by conducting 10 000 pseudo-experiments, separately for the 
muon and electron channels, using simulated tt and background events. The average devia¬ 
tions between extracted mass and JSF and their input values are obtained as a function of mf en 
and the bias is fit with a linear function for each generated JSF value. From these fits, additional 
small corrections for calibrating the top quark mass and the jet energy scale are derived as lin¬ 
ear functions of both the extracted top quark mass and JSF. The corrections are approximately 
—0.2 GeV for m t and —0.4% for the JSF. The statistical uncertainties of the method are also cor¬ 
rected by factors of approximately 1.04 that are derived from the widths of the corresponding 
pull distributions. 

The 2D ideogram fit to the combined electron and muon channels yields: 

m\ D = 172.14 ± 0.19 (stat+JSF) GeV, 

JSF 2D = 1.005 ± 0.002 (stat). 

As mt and the JSF are measured simultaneously, the statistical uncertainty in m t combines the 
statistical uncertainty arising from both components of the measurement. The uncertainty of 
the measurement agrees with the expected precision obtained by performing pseudo-experiments. 

The results in the individual muon and electron channels are compatible within their statistical 
uncertainties: 

y+jets: m\ D = 172.03 ± 0.27 (stat+JSF) GeV, JSF 2D = 1.007 ± 0.003 (stat), 
e+jets: m 2D = 172.26 ± 0.28 (stat+JSF) GeV, JSF 2D = 1.003 ± 0.003 (stat). 

The ID and hybrid analyses give results of 

m t 1D = 172.56 ± 0.12 (stat) GeV, 
m t hyb = 172.35 ± 0.16 (stat+JSF) GeV, 
jSF hyb = 1.002 ± 0.001 (stat). 
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respectively. Both the 2D and hybrid results for the JSF (JSF 2D and JSF hyb ) are within 0.5% of 
one. The results for m t and the JSF are compared in Fig. [4j which shows the two-dimensional 
statistical likelihoods obtained from data in the 2D and hybrid cases and m t from the ID anal¬ 
ysis. 

5.2 The all-jets channel 

As in the lepton+jets channel, a kinematic fit Il47l is used to improve the resolution of the recon¬ 
structed quantities and to check the compatibility of an event with the tt hypothesis. For each 
event, the inputs to the fit are the four-momenta of the six leading jets. The fit constrains these 
to the hypothesis of the production of two heavy particles of equal masses, each one decaying 
to a W boson with its invariant mass constrained to 80.4 GeV HTI and a bottom quark. The two 
b-tagged jets are candidates for the bottom quarks in the tt hypothesis, while the four untagged 
jets serve as candidates for the light quarks of the W boson decays. This leads to six possible 
parton-jet assignments per event and the assignment that fits best to the tt hypothesis based 
on the x 1 of the kinematic fit is chosen. As final selection criteria, we require P go f >0.1 for the 
kinematic fit with three degrees of freedom, and the two b quark jets be separated in t]-(p space 
by AR bb > 2.0. These requirements select 7049 events for the mass measurement in data and 
the fraction of signal events f s \g increases from 14% to 61% based on the simulation. 

For simulated tt events, the parton-jet assignments are classified as correct permutations and 
wrong permutations. The correct permutation class is defined in the same way as for the lep¬ 
ton+jets channel (Section |5.1| l. The wrong permutations class consists of permutations where at 
least one quark from the tt decay is not unambiguously matched with a distance of A R < 0.3 
to any of the six selected jets. For correct permutations, which compose 42% of the selected tt 
events, the kinematic fit improves the resolution of the fitted values of m t from 13.8 to 7.5 GeV. 

The multijet background from QCD is modeled using a control sample that is obtained from 
data with the same event selection but without the b tagging requirement. While this sample 
has a small contamination of a few percent coming from signal events, these have no influence 
on the results for the background model. For each event, the kinematic selection is applied to 
all possible assignments of the six jets to the six quarks from the tt hypothesis. The best fitting 
assignment is chosen and the event is used to model the background if it fulfills the ^gof and 
AR b b criteria. The modeled background is compared to the background predicted by an event 
mixing technique |[48l . Both predictions are found to agree within their uncertainties that are 
derived from the validation of the methods on simulated multijet events. 

Figure [5] compares data to the expectation from simulated tt signal and background estimate 
from the data for m[ ]t , m^ co , P go f, and AP bb . 

The ideogram method is calibrated for each combination of the 7iif en and JSF values by con¬ 
ducting 10 000 pseudo-experiments. The average deviations between extracted mass and JSF 
and their input values are obtained as a function of mf" and the bias is fit with a linear function 
for each generated JSF value. From these fits, additional small corrections for calibrating the 
top quark mass and the jet energy scale are derived as linear functions of both the extracted 
top quark mass and JSF. The corrections are approximately —0.6 GeV for m t and +1.0% for the 
JSF. The statistical uncertainties of the method are corrected by factors of approximately 1.005 
using values derived from the widths of the corresponding pull distributions. 

Applying the ideogram method on data with no prior knowledge on the JSF (2D), yields the 
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Figure 5: Distributions of (upper left) the reconstructed top quark mass from the kinematic fit, 
(upper right) the average reconstructed Wboson mass, (lower left) the goodness-of-fit probabil¬ 
ity, and (loser right) the separation of the two b quark jets for the all-jets channel. The simulated 
tt signal and the background from the control sample are normalized to data. The value of mo¬ 
used in the simulation is 172.5 GeV and the nominal jet energy scale is applied. The vertical 
bars show the statistical uncertainty and the hatched bands show the statistical and system¬ 
atic uncertainties added in quadrature. The lower portion of each panel shows the ratio of the 
yields between the collision data and the simulation. 
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CMS All-jets, 18.2 fb' 1 (8 TeV) 



Figure 6: The two-dimensional likelihood (—2Alog(£)) for the all-jets channel for the 2D, 
hybrid, and ID fits. The thick (thin) ellipses correspond to contours of —2Alog (£) = 1 (4) 
allowing the construction of the one (two) u statistical intervals of m t . For the ID fit, the thick 
and thin lines correspond to the one and two cr statistical uncertainties, respectively. 

results: 


m 2D = 17164 ± 0 32 ( stat+ j SF ) GeV/ 

JSF 2D = 1.011 ± 0.003 (stat). 

As m t and the JSF are measured simultaneously, the statistical uncertainty in m t combines the 
statistical uncertainty arising from both components of the measurement. The two additional 
free parameters in the fit, the signal fraction f s i g and the fraction of correct permutations / cp , 
are in agreement with the expectation from simulation. 

Using the JEC determined from y/Z+jet events in combination with the JSF prior from the 2D 
fit yields the results in the ID and hybrid approaches of 

m t 1D = 172.46 ± 0.23 (stat) GeV, 
m t hyb = 172.32 ± 0.25 (stat+JSF) GeV, 

JSF byb = 1.002 ± 0.001 (stat). 


For the all-jets channel, the fitted results for the JSF (JSF 2D and JSF hyb ) are both within 1.1% of 
one. While the JSF results from the 2D analyses in the lepton+jets and all-jets channels differ 
by 0.6%, the results from the hybrid analyses agree to within 0.2%. The all-jets results for m t 
and the JSF are compared in Fig. [6] which shows the two-dimensional statistical likelihoods 
obtained from data in the 2D and hybrid cases and m t from the ID analysis. 

5.3 The dilepton channel 

Figure[7]shows the distribution of mf MWT in data compared to a simulation with m t = 172.5 GeV 
for events containing at least one b jet. This channel is very clean with a negligible background 
from non tt sources and the collision and simulated events are in good agreement. 
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Figure 7: Distribution of »z( XMWT for the collision and simulated data with m t = 172.5 GeV. The 
vertical bars show the statistical uncertainty and the hatched bands show the statistical and 
systematic uncertainties added in quadrature. The lower section of the plot shows the ratio of 
the yields between the collision data and the simulation. 

AMWT masses are computed for all events in both the data and the simulations. The »z^ MWT 
distributions computed for each of the seven simulated tt mass samples are added to the distri¬ 
butions from the background samples, and these are treated as templates in a binned likelihood 
fit. To minimize the effects of any bias from the poorly populated tails of the distribution, we 
only examine events with m' XMW I between 100 and 400 GeV. For each of the seven mass tem¬ 
plates, a maximum likelihood fit is performed to the data distribution. A parabola is fit to the 
negative logarithms of the maximum likelihoods returned by the fits, and the minimum of the 
parabola is taken as the measured mass value. 

The fit is calibrated to correct for any biases induced by the reconstruction using pseudo-data. 
The calibration is performed by means of a test using the simulated templates for the top quark 
masses between 169.5 and 175.5 GeV. We randomly draw 1000 samples of events, each selected 
such that the total number of events is the same as in the full data sample. For each template, 
the 1000 measured masses are averaged together and subtracted from the input mass to obtain 
a numerical value for the bias induced by the fit. The bias is then parametrized as a linear 
dependence on the generated value of »z t , and the resulting calibration curve is used to correct 
for biases in the final result. 

The likelihoods obtained from a fit of each of the seven simulated templates to data, as well as 
a second-order polynomial fit to these values, are shown in Fig. [8] This yields an uncalibrated 
measured mass of m t = 172.77 ± 0.19 (stat) GeV. After correcting for the fit bias, the result for 
the top quark mass is found to be m t = 172.82 ± 0.19 (stat) GeV. 

The analysis was optimized with the value of ?zz t blinded. The optimization of the event se- 
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CMS Dilepton, 19.7 fb' 1 (8 TeV) 



m t [GeV] 


Figure 8 : Plot of the negative log-likelihood for data for the dilepton analysis. The continuous 
line represents a parabolic fit to the points. 

lection was done by minimizing the total expected (statistical+systematic) uncertainty This 
resulted in the restriction of the analysis to events containing only two b jets, rather than the 
requirement of at least one b jet which was used initially. 


6 Systematic uncertainties 

The systematic uncertainties affecting each of the measurements can be grouped into four dis¬ 
tinct categories: one experimental category and three theoretical categories that describe the 
modeling uncertainties. The experimental classification covers the uncertainties that arise from 
the precision of the calibration and resolution of the CMS detector and the effects coming from 
the backgrounds and pileup. The other three categories cover the modeling of the hard scatter¬ 
ing process and the associated radiation; non-perturbative QCD effects, such as the simulation 
of the underlying event and color reconnection; and the modeling of the light- and b-quark 
hadronization. Each of these is broken down into sub-categories leading to a total of 24 distinct 
systematic uncertainties. In each case the uncertainty is evaluated in terms of the largest shift 
that is observed in the value of m t that occurs when the parameter is varied by ±1 a, where a 
is the uncertainty assigned to that quantity. The only exception to this is if the statistical un¬ 
certainty in the observed shift is larger than the value of the calculated shift. In this case the 
statistical uncertainty is taken as the best estimate of the uncertainty in the parameter. 

6.1 Experimental effects 

• Intercalibration jet energy correction: This is the part of the JES uncertainty originat¬ 
ing from modeling of the radiation in the relative (pj- and (/-dependent) intercali¬ 
bration procedure. 

• In situ jet energy calibration: This is the part of the JES uncertainty coming from the 
uncertainties affecting the absolute JES determination using 7 /Z+jets events. 

• Uncorrelated jet energy correction: This is the uncertainty source coming from the 
statistical uncertainty in the in situ jet energy calibration, the contributions stemming 
from the jet energy correction due to pileup effects, the uncertainties due to the vari- 
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ations in the calorimeter response versus time, and some detector specific effects. To 
give a clear indication of the contribution to the JES uncertainty coming from pileup, 
we have sub-divided this uncertainty into non-pileup and pileup contributions. 

• Lepton energy scale (LES): Analogous to the JES, the energy scale of the leptons may 
also induce a systematic bias. A typical variation of 0.6% is taken for electrons in the 
barrel region and 1.5% in the detector endcaps. For the muons, the uncertainty is 
negligible. 

• E™ 1SS scale: Measurement of the Ej' ss is affected by the variation in LES and JES and 
by the uncertainty in scale of the unclustered energy The unclustered energy scale 
is varied independently of LES and JES to obtain the E™ ISS uncertainty 

• Jet energy resolution: The systematic uncertainty associated with the JER in the sim¬ 
ulation is determined by increasing or decreasing the JER by la. 

• b tagging: The uncertainty in the b tagging efficiency and misidentification probabil¬ 
ity of non-b jets may lead to varying background and signal levels. This uncertainty 
is estimated by varying the b tagging discriminator requirements in simulations. In 
the lepton+jets analysis, for example, the changes in the CSVM discriminator leads 
to an uncertainty in the b tagging efficiency of 1.2% and the false tagging rate of 
15%, both of which correspond to a la variation in the value of the b tagging scale 
factor. The all-jets and dilepton uncertainties are computed in a similar manner for 
the CSVT and CSVL discriminators, respectively Propagating the tagging efficiency 
uncertainty to the values of m t leads to the systematic uncertainty 

• Trigger: This systematic uncertainty captures the uncertainties related to the model¬ 
ing of the trigger efficiency, and is only significant for the all-jets measurement. 

• Pileup: During the data taking period, the instantaneous luminosity increased dra¬ 
matically during the year, leading to an increased number of simultaneous proton- 
proton interactions overlapping with the primary hard scattering (in-time pileup) 
and possible effects due to the detector response to previous collision events (out- 
of-time pileup). These effects are evaluated by using pseudo-experiments in which 
the average number of pileup events was varied by ±5%. 

• Backgrounds: The background contamination expected from simulation is <5% in 
the lepton+jets and dilepton channels. The effect of the background modeling on 
7?z t is estimated by varying the shape and normalization for each background within 
their uncertainty Uncertainties from simulated backgrounds are taken to be corre¬ 
lated across all the measurements. The only channel for which there is a significant 
non-tt background in the final fit sample is the all-jets channel. For this, the shape 
of the QCD multijet background is estimated from a control sample in the data. The 
method is validated using simulated QCD multijet events and with an alternative 
approach using event mixing in the data. The predicted background shapes are var¬ 
ied to cover the residual differences found in the validation. The uncertainties from 
the background estimation from control samples in the data are assumed to be un¬ 
correlated. 

• Fit calibration: For the calibration of the fits, the simulated samples are statistically 
limited. The uncertainty quoted is the statistical uncertainty in the residual bias in 
the fit calibration. 
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6.2 Theoretical and modeling uncertainties 

6.2.1 Hard scattering and radiation 

• Parton distribution functions: PDFs are used in modeling the hard scattering in 
proton-proton collisions in the simulations. The uncertainties in the PDFs and their 
effect on the measured value of m t are studied by reweighting a tt sample with 
different PDF eigenvectors using the PDF4LHC prescription Il49l . The reweighted 
events are used to generate pseudo-experiments and the variation in the fitted mass 
is quoted as the uncertainty l37l . 

• Renormalization and factorization scales: This uncertainty is estimated using the 
simultaneous variation of the renormalization and factorization scales by factors of 
2 and 0.5 in the matrix element calculation and the initial-state parton shower of the 
signal and the W+jets background. 

• ME-PS matching threshold: In the tt simulation, the matching thresholds used for 
interfacing the matrix elements (ME) generated with MadGraph to the PYTHIA 
parton showers (PS) are varied from the default of 40 GeV down to 30 GeV and up to 
60 GeV and the uncertainty is taken as the maximal difference in m t induced by this 
variation. 

• ME generator: The sensitivity to the parton-level modeling is estimated by com¬ 
paring the reference samples (MadGraph and pythia) to samples produced using 
POWHEG and pythia. The difference between the values of m t obtained with the 
two samples is taken as the uncertainty. 

• Top quark pj uncertainty: This term represents the uncertainty coming from the 
modeling of the top quark pj distribution in the ME generator. The uncertainty 
is estimated by taking the difference in shape between the parton level pj spectrum 
from the ME generator and the unfolded pj spectrum determined from the data f39l . 
The uncertainty is considered fully correlated across the measurements. 

6.2.2 Non-perturbative QCD 

• Underlying event: This represents the uncertainty in modeling the soft underlying 
hadronic activity in the event, which affects the simulation of both signal and back¬ 
ground. The uncertainty is estimated by comparing PYTHIA tunes with increased 
and decreased underlying event activity relative to a central tune. For this we com¬ 
pare the results for the Perugia 2011 tune to the results obtained using the Perugia 
2011 mpiHi and the Perugia 2011 Tevatron tunes Il50l . 

• Color reconnection: The effects of possible mismodeling of color reconnection are 
estimated by comparing the mass calculated using underlying event tunes with and 
without the inclusion of these effects. For these simulations the Perugia 2011 and 
Perugia 2011 no CR tunes are used Ii50l . The uncertainty is taken as the difference 
between the two computed values of nit- 

6.2.3 Hadronization 

• Flavor-dependent hadronization uncertainty: This is the part of the JES uncertainty 
that comes from differences in the energy response for different jet flavors and flavor 
mixtures with respect to those used in the calibration procedures. Four uncertain¬ 
ties are quoted that correspond to the uncertainties for light quarks (u, d, s), charm 
quarks, bottom quarks and gluons. These are evaluated by comparing Lund string 



19 


fragmentation (PYTHIA 6 Ifl6l ) and cluster fragmentation (HERWIG++ fl5ll ) for each 
category of jets. The models in pythia and herwig allow for the differences be¬ 
tween the jet types, and the uncertainty is determined by varying the jet energies 
within their respective flavor-dependent uncertainties. The full flavor-dependent 
uncertainty is obtained by taking a signed linear sum of these four contributions. 
For this we perform ± 1(7 shifts for each of the contributions and compute the total 
uncertainty from the sum of the ±1 a and — la shifts separately. As these are sym¬ 
metric, we quote the ± 1(7 shifts for the values of the uncertainties in Tables[ljj3] 

• b quark fragmentation and b hadron branching fraction uncertainties: This term pro¬ 
vides a description of the residual uncertainties not covered by the flavor-dependent 
hadronization term. It has two components: the uncertainty in the modeling of the 
b quark fragmentation function and the uncertainty from the measured b hadron 
semileptonic branching fractions. The b quark fragmentation function in PYTHIA 
is modeled using a Bowler-Lund model for the fragmentation into b hadrons. The 
fragmentation uncertainty is determined from the difference between a version tuned 
to ALEPH ||52| and DELPHI |53ll data and the pythia Z2* tune. Lastly, the uncer¬ 
tainty from the semileptonic b hadron branching fraction is obtained by varying by 
—0.45% and ±0.77%, which is the range of the measurements from B°/B 1 decays 
and their uncertainties HD . 


7 Individual channel results 

7.1 The lepton+jets channel 

After estimating the systematic uncertainties for the lepton+jets channel, the measurement of 
nit and the JSF from the 2D analysis gives 

m\ D = 172.14 ± 0.19 (stat+JSF) ± 0.59 (syst) GeV, 

JSF 2D = 1.005 ± 0.002 (stat) ± 0.007 (syst). 

The overall uncertainty in m t is 0.62 GeV and the measured JSF is compatible with the one 
obtained from events with Z bosons and photons |37| within the systematic uncertainties. 

The measurements from the ID and hybrid analyses are 

m} D = 172.56 ± 0.12 (stat) ± 0.62 (syst) GeV, 

m t hyb = 172.35 ± 0.16 (stat+JSF) ± 0.48 (syst) GeV. 

Thus the hybrid approach delivers the most precise measurement of the methods studied for 
the lepton+jets channel with a total uncertainty of 0.51 GeV. 

The breakdown of the systematic uncertainties for the three fits is shown in Table [l] In the 
lepton+jets and all-jets measurements several uncertainty sources yield opposite signs in the 
ID and 2D approaches. This arises because the untagged jets used for m(y C0 have a softer pj 
spectrum and larger gluon contamination compared to the b jets. As a consequence, the mea¬ 
surement of the JSF in the 2D measurement is more sensitive to low-px effects and radiation 
uncertainties than the ID measurement where the light-jet energies are bound to fulfill the 
W mass constraint. The net effect, when using a flat JSF, is that the uncertainties can be over¬ 
corrected in the 2D fit and thus their signs reverse. The hybrid fit makes optimal use of the 
available information and leads to partial cancelation of these uncertainties, resulting in the 
observed improvement of the precision of the mass measurement. 
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7 Individual channel results 


Table 1: Category breakdown of the systematic uncertainties for the 2D, ID, and hybrid mea¬ 
surements in the lepton+jets channel. Each term has been estimated using the procedures de¬ 
scribed in Section [6] The uncertainties are expressed in GeV and the signs are taken from the 
+1 cr shift in the value of the quantity Thus a positive sign indicates an increase in the value of 
7ni or the JSF and a negative sign indicates a decrease. For uncertainties determined on inde¬ 
pendent simulated samples the statistical precision of the shift is displayed. With the exception 
of the flavor-dependent JEC terms (see Section [6), the total systematic uncertainty is obtained 
from the sum in quadrature of the individual systematic uncertainties. 

Lepton+jets channel 


m t fit type 


Sm 2 t D (GeV) 

2D 

+JSF 

ID 

Sml D (GeV) 

hybrid 
6m h t yh (GeV) 

Experimental uncertainties 

Method calibration 

0.04 

0.001 

0.04 

0.04 

Jet energy corrections 





- JEC: Intercalibration 

<0.01 

<0.001 

+0.02 

+0.01 

- JEC: In situ calibration 

-0.01 

+0.003 

+0.24 

+0.12 

- JEC: Uncorrelated non-pileup 

+0.09 

-0.004 

-0.26 

-0.10 

- JEC: Uncorrelated pileup 

+0.06 

-0.002 

-0.11 

-0.04 

Lepton energy scale 

+0.01 

<0.001 

+0.01 

+0.01 

E™ 1SS scale 

+0.04 

<0.001 

+0.03 

+0.04 

Jet energy resolution 

-0.11 

+0.002 

+0.05 

-0.03 

b tagging 

+0.06 

< 0.001 

+0.04 

+0.06 

Pileup 

-0.12 

+0.002 

+0.05 

-0.04 

Backgrounds 

+0.05 

< 0.001 

+0.01 

+0.03 

Modeling of hadronization 

JEC: Flavor-dependent 





- light quarks (u d s) 

+0.11 

-0.002 

-0.02 

+0.05 

- charm 

+0.03 

<0.001 

-0.01 

+0.01 

- bottom 

-0.32 

<0.001 

-0.31 

-0.32 

- gluon 

-0.22 

+0.003 

+0.05 

-0.08 

b jet modeling 





- b fragmentation 

+0.06 

-0.001 

-0.06 

<0.01 

- Semileptonic b hadron decays 

-0.16 

<0.001 

-0.15 

-0.16 

Modeling of perturbative QCD 

PDF 

0.09 

0.001 

0.06 

0.04 

Ren. and fact, scales 

+0.17 + 0.08 

-0.004 ± 0.001 

-0.24 + 0.06 

-0.09 ± 0.07 

ME-PS matching threshold 

+0.11+0.09 

-0.002 + 0.001 

-0.07 + 0.06 

+0.03 ± 0.07 

ME generator 

-0.07 + 0.11 

-0.001 +0.001 

-0.16 + 0.07 

-0.12 + 0.08 

Top quark p-\ 

+0.16 

-0.003 

-0.11 

+0.02 

Modeling of soft QCD 

Underlying event 

+0.15 + 0.15 

-0.002 + 0.001 

+0.07 + 0.09 

+0.08 + 0.11 

Color reconnection modeling 

+0.11+0.13 

-0.002 + 0.001 

-0.09 ± 0.08 

+0.01 ± 0.09 

Total systematic 

0.59 

0.007 

0.62 

0.48 

Statistical 

0.20 

0.002 

0.12 

0.16 

Total 

0.62 

0.007 

0.63 

0.51 
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7.2 The all-jets channel 

The 2D analysis in the all-jets channel yields a measurement of 

ml D = 171.64 ± 0.32 (stat+JSF) ± 0.95 (syst) GeV, 

JSF 2D = 1.011 ± 0.003 (stat) ± 0.011 (syst), 

giving an overall uncertainty in the mass of 1.00 GeV. 

The measurements from the ID and hybrid analyses are 

m t 1D = 172.46 ± 0.23 (stat) ± 0.62 (syst) GeV, 
m t hyb = 172.32 ± 0.25 (stat+JSF) ± 0.59 (syst) GeV, 
with overall uncertainties of 0.66 and 0.64 GeV for the ID and hybrid fits, respectively 
The breakdown of the systematic uncertainties for the three fits is shown in Table [2] 

7.3 The dilepton channel 

For the dilepton channel the systematic uncertainties are defined as the difference between 
measurements of m t from pseudo-data events, selected at random from the MC events in the 
m t = 172.5 GeV template. For each category of systematic uncertainty, modified templates were 
produced with a given systematic variable shifted, generically by ±1 a. The fit is repeated using 
the modified pseudo-data and the respective mean is subtracted from the mean of the default 
tt MC simulation to calculate the final systematic uncertainty for each category This yields a 
final mass measurement of 


m t = 172.82 ±0.19 (stat) ±1.22 (syst) GeV. 

The breakdown of the systematic uncertainty for the dilepton mass measurement is shown in 
Table [3j In comparison with the lepton+jets (Table [l]l and the all-jets (Table [2j channels, the 
systematic uncertainties are similar in size with the exception of the factorization and renor¬ 
malization, and b fragmentation terms, both of which are significantly larger. Studies of these 
indicate that this is probably the result of an increased boost of the visible decay products, 
coupled to the weak constraint of the E™ 1SS on the energies of the two neutrinos. 

7.4 The 2010 and 2011 measurements 

The published CMS measurements are based on i/s = 7TeV data recorded during 2010 and 
2011. Although much less precise than the new measurements, they come from independent 
data sets and have different sensitivities to the various systematic uncertainties. These are 
included in the combined mass analysis, which is discussed in Section [9] For completeness 
we summarize these measurements in Table [4] below. The analysis techniques used for each 
of these are very similar to those used for the 2012 analyses. The dilepton results both use the 
AMWT method, which is described in Section |4.2[ and the lepton+jets (all-jets) result comes 
from the 2D (ID) ideogram technique, which is described in Section [4~Tj 

8 Measured top quark mass as a function of kinematic observ¬ 
ables 

To search for possible biases in our measurements and the potential limitations of current event 
generators, a series of differential measurements of m t as a function of the kinematic properties 
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Table 2: Category breakdown of the systematic uncertainties for the 2D, ID and hybrid mea¬ 
surements in the all-jets channel. Each term has been estimated using the procedures described 
in Section[6j The uncertainties are expressed in GeV and the signs are taken from the +1(7 shift 
in the value of the quantity. Thus a positive sign indicates an increase in the value of m t or 
the JSF and a negative sign indicates a decrease. For uncertainties determined on independent 
simulated samples the statistical precision of the shift is displayed. With the exception of the 
flavor-dependent JEC terms (see Section [6), the total systematic uncertainty is obtained from 
the sum in quadrature of the individual systematic uncertainties. 

All-jets channel 


rn t fit type 


(5m t 2D (GeV) 

2D 

()JSF 

ID 

Sm t 1D (GeV) 

hybrid 
6m byb , (GeV) 

Experimental uncertainties 

Method Calibration 

0.06 

0.001 

0.06 

0.06 

Jet energy corrections 





- JEC: Intercalibration 

<0.01 

<0.001 

+0.02 

+0.02 

- JEC: In situ calibration 

-0.01 

<0.001 

+0.23 

+0.19 

- JEC: Uncorrelated non-pileup 

+0.06 

-0.001 

-0.19 

-0.16 

- JEC: Uncorrelated pileup 

+0.04 

<0.001 

-0.08 

-0.06 

Jet energy resolution 

-0.10 

+0.001 

+0.03 

+0.02 

b tagging 

+0.02 

<0.001 

+0.01 

+0.02 

Pileup 

-0.09 

+0.002 

+0.02 

<0.01 

Backgrounds 

-0.61 

-0.007 

-0.14 

-0.20 

Trigger 

+0.04 

<0.001 

-0.01 

<0.01 

Modeling of hadronization 

JEC: Flavor-dependent 





- light quarks (u d s) 

+0.10 

-0.001 

-0.02 

+0.00 

- charm 

+0.03 

-0.001 

-0.01 

-0.01 

- bottom 

-0.30 

+0.000 

-0.29 

-0.29 

- gluon 

-0.17 

+0.002 

+0.02 

-0.02 

b jet modeling 





- b fragmentation 

+0.08 

-0.001 

+0.03 

+0.04 

- Semileptonic b hadron decays 

-0.14 

<0.001 

-0.13 

-0.13 

Modeling of perturbative QCD 

PDF 

0.06 

<0.001 

0.03 

0.03 

Ren. and fact, scales 

+0.29 + 0.16 

-0.005 ± 0.001 

-0.19 + 0.11 

-0.12 + 0.12 

ME-PS matching threshold 

+0.18 + 0.16 

-0.002 ± 0.001 

+0.12 + 0.11 

+0.13 + 0.12 

ME generator 

-0.04 ± 0.20 

-0.002 ± 0.002 

-0.18 + 0.14 

-0.16 + 0.14 

Top quark p-\ 

+0.04 

+0.001 

+0.08 

+0.06 

Modeling of soft QCD 

Underlying event 

+0.27 ± 0.25 

-0.002 ± 0.002 

+0.13 + 0.18 

+0.14 + 0.18 

Color reconnection modeling 

+0.35 ± 0.22 

-0.003 ± 0.002 

+0.14 + 0.16 

+0.16 + 0.16 

Total systematic 

0.95 

0.011 

0.62 

0.59 

Statistical 

0.32 

0.003 

0.23 

0.25 

Total 

1.00 

0.011 

0.66 

0.64 


of the ft system is performed. To maximize the accuracy of the results, the study is performed 
in the lepton+jets channel using the hybrid fit technique. The variables are chosen to probe 
potential effects from color reconnection, initial- and final-state radiation, and the kinematics 
of the jets coming from the top quark decays. 
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Table 3: Category breakdown of the systematic uncertainties for the AMWT measurement in 


the dilepton channel. Each term has been estimated using the procedures described in Sec¬ 
tion^ The uncertainties are expressed in GeV and the signs are taken from the +1(7 shift in the 
value of the quantity. Thus a positive sign indicates an increase in the value of /n t and a nega¬ 
tive sign indicates a decrease. For uncertainties determined on independent simulated samples 
the statistical precision of the shift is displayed. With the exception of the flavor-dependent JEC 
terms (see Section [6j, the total systematic uncertainty is obtained from the sum in quadrature 
of the individual systematic uncertainties. 

Dilepton channel 

Sm t (GeV) 

Experimental uncertainties 

Method calibration 

0.03 

Jet energy corrections 


- JEC: Intercalibration 

+0.03 

- JEC: In situ calibration 

+0.24 

- JEC: Uncorrelated non-pileup 

-0.28 

- JEC: Uncorrelated pileup 

-0.12 

Lepton energy scale 

+0.12 

Ej' ss scale 

+0.06 

Jet energy resolution 

+0.06 

b tagging 

+0.04 

Pileup 

+0.04 

Backgrounds 

+0.02 

Modeling of hadronization 

JEC: Flavor-dependent 


- light quarks (u d s) 

+0.02 

- charm 

+0.02 

- bottom 

-0.34 

- gluon 

+0.06 

b jet modeling 


- b fragmentation 

-0.69 

- Semileptonic b hadron decays 

-0.17 

Modeling of perturbative QCD 

PDF 

0.16 

Ren. and fact, scales 

-0.75 ± 0.20 

ME-PS matching threshold 

-0.12 + 0.20 

ME generator 

-0.24 ± 0.20 

Top quark pj 

-0.25 

Modeling of soft QCD 

Underlying event 

+0.04 ± 0.20 

Color reconnection modeling 

-0.11 + 0.20 

Total systematic 

1.22 

Statistical 

0.19 

Total 

1.24 


For each measurement, the hybrid analysis method is applied to subsets of events defined 
according to the value of a given kinematic event observable after the kinematic fit. The con¬ 
tribution of the external JSF constraint is fixed to 50% to ensure consistency of all bins with 
the inclusive result. Constant shifts in the measured m t values may arise due to the systematic 
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8 Measured top quark mass as a function of kinematic observables 


Table 4: CMS measurements of the top quark mass using the data recorded at a/s = 7 TeV. 



Analysis 

Reference 

mt 

(GeV) 

Stat. uncertainty 
(GeV) 

Syst. uncertainty 
(GeV) 

2010 

dilepton (AMWT) 

El 

175.50 

4.60 

4.52 

2011 

lepton+jets (2D) 

m 

173.49 

0.27 

1.03 

2011 

all-jets (ID) 

m 

173.49 

0.69 

1.23 

2011 

dilepton (AMWT) 

m 

172.50 

0.43 

1.46 


uncertainties of the inclusive measurement or from the use of different ni t values in data and 
simulations. To search for kinematics-dependent biases the value of the mean measured top 
quark mass is subtracted and the results are expressed in the form m t — (nit), where the mean 
comes from the inclusive measurement on the specific sample. In each case, the event sample 
is divided into 3 to 5 bins as a function of the value of the kinematic observable and we popu¬ 
late each bin using all permutations which lie within the bin boundaries. As some observables 
depend on the jet-quark assignment that cannot be resolved unambiguously, such as the pj of 
a reconstructed top quark, a single event is allowed to contribute to multiple bins. 

To aid in the interpretation of a difference between the value of nit ~ (nit) arid the prediction 
from a simulation in the same bin, a bin-by-bin calibration of the results is performed using 
the MadGraph + PYTHIA simulation. This is performed using the same technique as for the 
inclusive measurement Il55l except that it is performed on each bin separately. Thus, after 
calibration the value in each bin can be interpreted in terms of its agreement with respect to the 
inclusive measurement. 


For eight kinematic variables the results for the calibrated mass difference, m t — (nit), are shown 
as a function of the chosen variable, and we compare the results to the predictions of seven 
different simulations. For each plotted point the statistical uncertainty and the dominant sys¬ 
tematic uncertainties are combined in quadrature, where the latter include the JES (pj-, rj- and 
flavor-dependent), JER, pileup, b fragmentation, renormalization and factorization scales, and 
ME-PS matching threshold. The systematic uncertainties are assumed to be correlated among 
all bins, so that any constant shift is removed by subtracting (nit). We note that this approxi¬ 
mation may underestimate the uncertainties from the pj ///-dependent JES. 


For each plot we compare the data to simulations based on LO (MadGraph and SHERPA) 
and NLO (POWHEG and MC@NLO) matrix element calculations with both string (PYTHIA) and 
cluster (HERWIG) models for fragmentation. We also vary the choice of underlying event tune 
from Z2* to Perugia 2011 both with and without color reconnections, and the AUET2 tune. With 
the exception of the MC@NLO and SHERPA simulations, which are only used for this study, these 
are the same simulation as those discussed in Section 3.1 The simulations used for this study 
are 


• MadGraph with the PYTHIA Z2* tune, which is the simulation used in the mass 
determinations lH4HT6ll56ll ; 

• MadGraph with the PYTHIA Perugia 2011 tune lfl4l|T6ll50B ; 

• MadGraph with the PYTHIA Perugia 2011 noCR tune ITTlTfi , 501; 

• SHERPA 1.4.0 with up to 4 additional jets from the LO matrix element Il57ll58ll : 

• POWHEG with the PYTHIA Z2* tune lH6ll2Tti25ll56l ; 

• POWHEG with the HERWIG 6.520 AUET2 tune I2lti25ll59l ; 

• MC@NLO 3.41 with the HERWIG 6.520 default tune l59l - ?6ll . 
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Table 5: Comparison of different simulations and the data. The summed x 2 values and number 
of standard deviations are computed for the 27 points entering Figs. [9] and 10 assuming two- 
sided Gaussian statistics. 


Simulation 

x z 

Standard deviations 

MG + PYTHIA 6 Z2* 

17.55 

0.10 

MG + PYTHIA 6 Pll 

37.68 

1.73 

MG + PYTHIA 6 PllnoCR 

31.57 

1.15 

POWHEG + PYTHIA 6 Z2* 

19.70 

0.20 

POWHEG + HERWIG 6 

76.48 

4.84 

MC@NLO + HERWIG 6 

20.47 

0.24 

SHERPA 

46.79 

2.56 


The variables were chosen for their potential sensitivity to modeling the kinematics of top 
quark production (Fig. [9jl and decay (Fig. [To| . No significant deviation in the value of the 
measured m t is observed, indicating that within the current precision, there is no evidence for 
a bias in the measurements. The agreement between the data and each of the simulations is 
quantified in Table [ 5 ] Here we show the cumulative x 2 f° r the 27 degrees of freedom repre¬ 
sented by the eight distributions studied (Figs. [9] and [To]) and the corresponding number of 
standard deviations between the data and the simulation, where we have assumed two-sided 
Gaussian confidence intervals for each simulation. In all cases, with the possible exception of 
POWHEG + HERWIG 6 simulation, the data is well described by the models. 


9 Combining the mass measurements 

In this section, results for the combined top quark mass measurement are presented. As in¬ 
puts we use the new results presented in this paper and the published CMS measurements 
from the 2010 Il54ll and 2011 Il43l l48l l55l analyses. To combine the results, the best linear unbi¬ 
ased estimate method (BLUE) Il62l is used. This determines a linear combination of the input 
measurements which takes into account statistical and systematic uncertainties by minimizing 
the total uncertainty of the combined result. The procedure takes account of the correlations 
that exist between the different uncertainty sources through the use of correlation coefficients. 
These are chosen to reflect the current knowledge of the uncertainties for both the correlations 
between measurements in a given decay channel from different years (p c han) and between the 
measurements in different decay channels from the same year (p ye ar)- The nominal values are 
set to either zero for uncorrelated or unity for fully correlated (see Table [6|. Because the mea¬ 
surements from the 2012 analyses are significantly more precise, both statistically and system¬ 
atically, than those from the 2010 and 2011 analyses, the use of unity coefficients for p c h an and 
Pyear is problematic. To mitigate this, we have chosen to perform combinations in which the 
correlation coefficients are limited to value of less than unity. This has been done by setting the 
correlation coefficients for each pair of measurements in the fully correlated cases to p = Jj/ CTj, 
where cy and cry are the uncorrelated components of the uncertainties in measurements i and /, 
respectively, and G < cry. For all of the measurements, the statistical uncertainties are assumed 
to be uncorrelated. 

9.1 Measurement permutations 

The precision of any combination of the measurements will be dominated by the set of new 
measurements, derived from the 2012 data. To investigate the effect of the choice of fit method 
on the result, we perform a series of combinations in which the 2012 inputs from the lepton+jets 
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9 Combining the mass measurements 
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Figure 9: Measurements of tih as a function of the transverse momentum of the hadronically 
decaying top quark (pj Wd ), the invariant mass of the tt system (m tt ), the transverse momentum 
of the tt system (pj), and the number of jets with pj > 30 GeV. The filled circles represent the 
data, and the other symbols are for the simulations. For reasons of clarity the horizontal error 
bars are shown only for the data points and each of the simulations is shown as a single offset 
point with a vertical error bar representing its statistical uncertainty. The statistical uncertainty 
of the data is displayed by the inner error bars. For the outer error bars, the systematic uncer¬ 
tainties are added in quadrature. The open circles correspond to MadGraph with the PYTHIA 
Z2* tune, the open squares to MadGraph with the pythia Perugia 2011 tune, and the open 
triangles represent MadGraph with the PYTHIA Perugia 2011 noCR tune. The open diamonds 
correspond to POWHEG with the PYTHIA Z2* tune and the open crosses correspond to POWHEG 
with HERWIG 6. The filled stars are for MC@NLO with HERWIG 6 and the open stars are for 
SHERPA. 


and all-jets decay channels are varied. For simplicity of discussion, these are classified accord¬ 
ing to the type of fit used for each channel. They are labeled as follows: 2 for a 2D fit, 1 for a ID 
or AMWT fit, and h for a hybrid fit. Thus a lepton+jets:all-jets:dilepton fit is denoted 211 in the 
case of a 2D fit for the lepton+jets channel a ID fit for the all-jets channel and the AMWT fit in 
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Figure 10: Measurements of m t as a function of the pj of the b jet assigned to the hadronic decay 
branch (p j had ), the pseudorapidity of the b jet assigned to the hadronic decay branch (| p b,had |), 
the AR between the b jets (AR bb ), and the AR between the light-quark jets (AR q q). The filled cir¬ 
cles represent the data, and the other symbols are for the simulations. For reasons of clarity the 
horizontal error bars are shown only for the data points and each of the simulations is shown 
as a single offset point with a vertical error bar representing its statistical uncertainty. The sta¬ 
tistical uncertainty of the data is displayed by the inner error bars. For the outer error bars, the 
systematic uncertainties are added in quadrature. The open circles correspond to MadGraph 
with the PYTHIA 72* tune, the open squares to MadGraph with the PYTHIA Perugia 2011 
tune, and the open triangles represent MadGraph with the PYTHIA Perugia 2011 noCR tune. 
The open diamonds correspond to POWHEG with the PYTHIA 72* tune and the open crosses 
correspond to POWHEG with HERWIG 6. The filled stars are for MC@NLO with HERWIG 6 and 
the open stars are for SHERPA. 


the dilepton channel. 

The most precise set of non-hybrid measurements corresponds to the set 211, which gives a 


































28 


9 Combining the mass measurements 


Table 6: Nominal correlation coefficients for the systematic uncertainties. The term p c h an is the 
correlation factor for measurements in the same top quark decay channel, but different years 
and the term p yeai is the correlation between measurements in different channels from the same 
year. 


Correlations 



Pchan 

Pyear 

Experimental uncertainties 

Method calibration 

0 

0 

JEC: Intercalibration 

1 

i 

JEC: In situ calibration 

1 

i 

JEC: Uncorrelated non-pileup 

0 

i 

Lepton energy scale 

1 

i 

E™ ss scale 

1 

i 

Jet energy resolution 

1 

i 

b tagging 

1 

i 

Pileup 

0 

i 

Non-tt background (data) 

0 

0 

Non-tt background (simulation) 

1 

1 

Trigger 

0 

0 

Modeling of hadronization 

JEC: Flavor-dependent 

1 

1 

b jet modeling 

1 

1 

Modeling of perturbative QCD 

PDF 

1 

1 

Ren. and fact, scales 

1 

1 

ME-PS matching threshold 

1 

1 

ME generator 

1 

1 

Top quark pj 

1 

1 

Modeling of soft QCD 

Underlying event 

1 

1 

Color reconnection modeling 

1 

1 


result of 

m t = 172.40 ± 0.13 (stat+JSF) ± 0.54 (syst) GeV (211 combination). 

To verify that this gives the most precise combination, combinations are performed using the 
other permutations of the 2012 measurements. The results, listed in Table [7j are in good agree¬ 
ment with the 211 result but have less precision, as expected. 

For the hybrid results, the effect of constraining the JSF factor in the mass fits can be examined. 
There are three significant new permutations to consider, the hll, hhl, and 2hl combinations. 
The results, shown in Table |7j are in good agreement with the 211 result, with the hll and hhl 
combinations giving the most precise measurements, as expected. For these the results are 

mt = 172.45 ± 0.13 (stat+JSF) ± 0.47 (syst) GeV (hll combination), 
nit = 172.44 ± 0.13 (stat+JSF) ± 0.47 (syst) GeV (hhl combination), 

both with an overall improvement in precision of 0.07 GeV with respect to the 211 analysis, and 
a total uncertainty of 0.48 GeV. 


















9.2 Anticorrelation effects 
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Table 7: Combination results for the permutations of the 2D, ID, and hybrid measurements. 
The permutation order is defined to be lepton+jets:all-jets:dilepton, thus 211 corresponds to the 
2D lepton+jets:lD all-jets:AMWT dilepton combination. 


Combination 

m t 

(GeV) 

Stat+JSF uncertainty 
(GeV) 

Syst uncertainty 
(GeV) 

211 

172.40 

0.13 

0.54 

121 

172.61 

0.11 

0.57 

221 

172.30 

0.15 

0.58 

111 

172.66 

0.12 

0.56 

hll 

172.45 

0.13 

0.47 

hhl 

172.44 

0.13 

0.47 

2 hl 

172.35 

0.14 

0.53 


9.2 Anticorrelation effects 


For the results presented here, the signs of most of the uncertainty contributions are well de¬ 
fined (i.e. for a 1<7 shift in a given quantity, the statistical component of the estimated systematic 
uncertainty is significantly smaller than the value of the uncertainty). This allows a compari¬ 
son of the signs of the systematic uncertainties for the different channels and for the different 
fitting techniques. An anticorrelation (i.e. opposite signs) is observed between several of the 
terms when comparing the results from a 2D and a ID (or AMWT) fit. However, if the 2D fit 
is replaced by the corresponding hybrid result, the anticorrelations are removed. This is illus¬ 
trated in Fig. Ill] which shows the uncertainty correlations between the lepton+jets and all-jets 


channels for the 2D vs. ID and the hybrid vs. ID cases. In the 2D vs. ID plot (Fig. 11 left) we 
observe a significant number of anticorrelated terms (coming primarily from the JES and pileup 
terms), whereas in the hybrid vs. ID plot (Fig. [Tl] right) we see no significant anti-correlations. 
Given the uncertainty terms that vary between the 2D and hybrid treatments, it is believed that 
the observed effect arises from the variation in the JSF factors between the 2D, ID, and hybrid 
results (Sections |5. 1 1 and px2[ . 


These effects are not considered in the standard 211 combination as the input correlation coeffi¬ 
cients are positive for all of the correlated cases (see Table [6]. To estimate the effect of including 
anticorrelations, the correlation coefficients are set to negative values for the cases where an 
anticorrelation (opposite sign) is observed and positive values where a positive (same sign) 
or neutral (statistically limited) correlation is observed and the 211 combination analysis is re¬ 
peated. This gives a result of 172.40 ± 0.13 (stat+JSF) ± 0.47 (syst)GeV. Thus, while the result 
for the mass is unchanged, the systematic uncertainty is decreased and becomes comparable to 
that achieved in the hybrid combinations. 


10 Results 

Based on the expected uncertainties for each of the individual measurements (Tables [l]-[3]l and 
the consistency of the hybrid and ID results for the JSF (Sections |5.1||5.2) , the hhl combination 
is chosen as the preferred result. Combining the seven input measurements (four from y/s = 
7 TeV and three from this analysis) gives a combined top quark mass measurement of 

m t = 172.44 ± 0.13 (stat+JSF) ± 0.47 (syst) GeV, 

for which the combination x 2 is 2.5 for six degrees of freedom, corresponding to a probability 
of 87%. This is compared to the full set of Run 1 measurements in Fig. [12] where the current 
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10 Results 


CMS 8 TeV CMS 8 TeV 



Figure 11: Systematic uncertainty correlations for mass measurements in the lepton+jets and 
all-jets channels. Each point represents a single systematic uncertainty taken from Tables [l] 
and [2j Panel (left) for the 2D lepton+jets and ID all-jets measurements, and Panel (right) for 
the hybrid lepton+jets and the ID all-jets measurements. The filled circles correspond to the 
systematic uncertainties which show a positive correlation between the two fit methods and the 
open circles to the systematic terms which show a negative correlation. The points shown as 
filled squares are those for which the systematic estimation is dominated by a statistical uncer¬ 
tainty, so no clear categorization is possible. The vertical and horizontal error bars correspond 
to the statistical uncertainties in the systematic uncertainties. 

world average 13 and Tevatron 0 combinations are also shown. The result is consistent with 
all of the published LHC measurements and is the most precise measurement to date with a 
precision of 0.3%. 


Table 8: Correlations between input measurements. The elements in the table are labelled ac¬ 
cording to the analysis they correspond to (rows and columns read as 2010,2011,2012 followed 
by the tt decay channel name). 




2010 

2011 

2012 

dilepton 

dilepton 

lepton+jets 

all-jets 

dilepton 

lepton+jets 

all-jets 

2010 

dilepton 

1.00 








dilepton 

0.15 

1.00 






2011 

lepton+jets 

0.09 

0.37 

1.00 






all-jets 

0.10 

0.62 

0.31 

1.00 





dilepton 

0.09 

0.26 

0.17 

0.17 

1.00 



2012 

lepton+jets 

0.05 

0.21 

0.30 

0.26 

0.26 

1.00 



all-jets 

0.06 

0.20 

0.27 

0.28 

0.32 

0.61 

1.00 


The correlations between each of the measurements is shown in Table |8j Figure 13 shows the 
combination coefficients and pulls, where the pull is defined as — m™ p as ) / v c^ eas — u? om 

where m {:°™ b and mjjjp 15 are the combined and the individual measurements of m t/ respectively, 
and c 7 tomb and <7 mea s are the corresponding total uncertainties. The 2010 measurement con¬ 
tributes very little to the overall result. As the treatment of the systematic uncertainty for this 
analysis is the least sophisticated of the seven measurements, the final combination is repeated 
to verify that it does not influence the final result. Excluding this measurement produces neg¬ 
ligible changes in the values of m t or its total uncertainty, Sm t . For the combination of the 
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1 1 1 1 1 1 


i | i i i i | i i 

CMS 2010, dilepton 


w 175.50 ±4.60 ±4.60 GeV 

JHEP 07 (2011) 049, 36 pb 1 


(value ± stat ± syst) 

CMS 2011, dilepton 


172.50 ± 0.43 ± 1.43 GeV 

EPJC 72 (2012) 2202, 5.0 fb 


(value ± stat ± syst) 

CMS 2011, all-jets 


w 173.49 ± 0.69 ± 1.21 GeV 

EPJC 74 (2014) 2758, 3.5 fb' 1 


(value ± stat ± syst) 

CMS 2011, lepton+jets 


9 173.49 ± 0.43 ± 0.98 GeV 

JHEP 12 (2012) 105, 5.0 fb >' 


(value ± stat ± syst) 

CMS 2012, dilepton 

■ 

172.82 ± 0.19 ± 1.22 GeV 

This analysis, 19.7 fb' 1 


(value ± stat ± syst) 

CMS 2012, all-jets 


172.32 ± 0.25 ± 0.59 GeV 

This analysis, 18.2 fb' 1 


(value ± stat ± syst) 

CMS 2012, lepton+jets 


172.35 ± 0.16 ± 0.48 GeV 

This analysis, 19.71b' 1 

\SS: 

(value ± stat ± syst) 

CMS combination 

♦ 

172.44 ±0.13 ±0.47 GeV 



(value ± stat ± syst) 

Tevatron combination (2014) 


. . 

arXiv: 1407.2682 


174.34 - 0.37 ± 0.52 GeV 



(value ± stat ± syst) 

World combination 2014 



ATLAS. CDF. CMS. DO 


173.34 ±0.27 ±0.71 GeV 

arXiv:1403.4427 

_i__i_i_i_i__i_ 

_ 

(value ± stat ± syst) 

_1__1_1_1_1__1_1_ 


165 170 175 180 

m t [GeV] 


Figure 12: Summary of the CMS m t measurements and their combination. The thick error bars 
show the statistical uncertainty and the thin error bars show the total uncertainty Also shown 
are the current Tevatron |8| and world average 0 combinations. 

remaining six measurements the y 2 is 2.3 for five degrees of freedom, corresponding to a prob¬ 
ability of 80%. 

The breakdown of the systematic uncertainties for the combination is shown in Table [9] The 
dominant uncertainty in the measurement arises from the modeling of the hadronization, with 
0.33 GeV coming from the flavor-dependent jet energy corrections and a further 0.14 GeV com¬ 
ing from the b jets. There are a further six terms with uncertainties in the range of 0.11-0.12 GeV. 
Of these, four are coming from theory and only two, the JEC in situ (0.12 GeV) and the JEC Un¬ 
correlated non-pileup (0.10 GeV) are experimental. The theoretical uncertainties are computed 
using the same models so they should be fully correlated. For the two experimental terms, 
the strength of the assumed correlations is varied by 50% of their nominal values to check the 
sensitivity to the assumed correlation strength. In both cases this produces changes of less than 
0.01 GeV in m t and 5m t . We therefore conclude that the result is quite stable against reasonable 
changes in the assumed correlation strength. 

Although we do not believe that the use of 100% correlation strengths is appropriate to use 
for the correlated systematic uncertainties, for completeness we have re-run the final combi¬ 
nation without the constraint on the correlation strengths. In this case we observe shifts of 
—0.28 GeV in m t and —0.03 GeV in 5m t . For this combination, four of the seven measurements 
have negative combination coefficients and the central mass lies outside of the boundaries of 
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CMS 


-2 0 2 

Pull 

Figure 13: Results of the BLUE combining procedure on the CMS measurements showing (left) 
the combination coefficients, and (right) the pulls for each contribution. 

the measurements. This corresponds to the result obtained using the standard BLUE method. 

Figure [14] shows the mass values obtained from each of the three channels separately. These 
correspond to combinations h2 (2012, 2011) for the lepton+jets channel. 111 (2012, 2011, 2010) 
for the dilepton channel, and hi (2012, 2011) for the all-jets channel, respectively. The results 
are all in good agreement with the combined measurement. 


1 

2010 dilepton 


1 

< 0.05 % 

2011 di lepton 


1.1 % 

2011 all-jets 


0.2 % 

2011 lepton+jets 


6.6 % 

2012 di lepton 


3.1 % 

2012 all-jets 

■ 

16 . 6 % 

2012 lepton+jets 

72.5 % 

_1 _ 

_ 1 _ 


-100 0 100 


BLUE Combination Coefficient [%] 


CMS 


T 


2010 dilepton 

2011 dilepton 
2011 all-jets 

2011 lepton+jets 

2012 dilepton 
2012 all-jets 
2012 lepton+jets 

_I_ 


0.48 

0.04 

0.80 

1.11 

0.33 

- 0.28 

- 0.63 


11 Summary 

A new set of measurements of the top quark mass has been presented, based on the data 
recorded by the CMS experiment at the LHC at i/s = 8TeV during 2012, and corresponding 
to a luminosity of 19.7 fb 1 . The top quark mass has been measured in the lepton+jets, all-jets 
and dilepton decay channels, giving values of 172.35 ± 0.16 (stat) ± 0.48 (syst)GeV, 172.32 ± 
0.25 (stat) ± 0.59 (syst) GeV, and 172.82 ± 0.19 (stat) ± 1.22 (syst) GeV, respectively. Individu¬ 
ally, these constitute the most precise measurements in each of the decay channels studied. 
When combined with the published CMS results at >/s = 7TeV, a top quark mass measure¬ 
ment of 172.44 ± 0.13 (stat) ± 0.47 (syst) GeV is obtained. This is the most precise measurement 
of m t to date, with a total uncertainty of 0.48 GeV, and it supersedes all of the previous CMS 
measurements of the top quark mass. 

The top quark mass has also been studied as a function of the event kinematical properties in 
the lepton+jets channel. No indications of a kinematical bias in the measurements is observed 
and the data are consistent with a range of predictions from current theoretical models of tt 
production. 
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Table 9: Category breakdown of systematic uncertainties for the combined mass result. The 
uncertainties are expressed in GeV. 


Combined m t result 

NjZt(GeV) 

Experimental uncertainties 

Method calibration 

0.03 

Jet energy corrections 


- JEC: Intercalibration 

0.01 

- JEC: In situ calibration 

0.12 

- JEC: Uncorrelated non-pileup 

0.10 

Lepton energy scale 

0.01 

Ef ss scale 

0.03 

Jet energy resolution 

0.03 

b tagging 

0.05 

Pileup 

0.06 

Backgrounds 

0.04 

Trigger 

<0.01 

Modeling of hadronization 

JEC: Flavor 

0.33 

b jet modeling 

0.14 

Modeling of perturbative QCD 

PDF 

0.04 

Ren. and fact, scales 

0.10 

ME-PS matching threshold 

0.08 

ME generator 

0.11 

Top quark pj 

0.02 

Modeling of soft QCD 

Underlying event 

0.11 

Color reconnection modeling 

0.10 

Total systematic 

0.47 

Statistical 

0.13 

Total Uncertainty 

0.48 
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